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Abstract We present a new solution for the rotation curves of galactic disks with 
gravitational potential of the Yukawa type. We follow the technique employed by 
Toomre in 1963 in the study of galactic disks in the Newtonian theory. This new 
solution allows an easy comparison between the Newtonian solution and the Yukawian 
one. Therefore, constraints on the parameters of theories of gravitation can be imposed, 
which in the weak field limit reduce to Yukawian potentials. We then apply our formulae 
to the study of rotation curves for a zero-thickness exponential disk and compare it with 
the Newtonian case studied by Freeman in 1970. As an application of the mathematical 
tool developed here, we show that in any theory of gravity with a massive graviton (this 
means a gravitational potential of the Yukawa type) , a strong limit can be imposed on 
the mass (m g ) of this particle. For example, in order to obtain a galactic disk with a 
scale length of b ~ lOkpc, we should have a massive graviton of m g << 10~ j9 g. This 
result is much more restrictive than those inferred from solar system observations. 

PACS 04.50. +h ■ 98.62.Hr 



1 Introduction 

General relativity agrees, as is well known, with Newtonian gravitation in the weak 
field limit (correspondence principle). Some alternatives to general relativity theory, 
however, not necessary respect this principle. In this way, independent of the specific 
details of the models or theories (for example, scalar-tensor theories of gravity, 
nonsymmetric gravitational theory, etc.), most of them converge to the same weak 
field limit: the gravitational potential is Yukawa-like (e.g., Refs. [l][2][3]|4]). 
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A first motivation of the present paper is to present a mathematical tool that can 
be used, independently of the particular theory of gravity, for the study of the rotation 
curves of galaxies. In this case, there is only one constrain: the potential must be 
Yukawa-like. 

Now, let us start recalling that the Newtonian potential 0, as is well known, follows 
the Poisson equation, namely 

V 2 = -4vrGp. (1) 



As a result, the potential of a point mass m at a distance r reads 

Gm 



(2) 



where G is the gravitational constant. 

Considering a gravitational potential of the Yukawa-type (hereafter named 
Yukawian gravitation potential), we have for a point mass m a potential given by 

Gm _ r / A , . 

4>= e , (3) 

r 

where A is a constant. Note that, for A — > oo this potential becomes identical to the 
Newtonian one. 

The field equation in this case reads 

V 2 - = -4nGp. (4) 

Recall that the constant A that appears in the Yukawa potential is the Compton 
wavelength of the exchange particle of mass m g . Using this interpretation for the 
present case one can think of the exchange particle is a massive graviton. In particular, 
the second motivation of the present paper is to verify if it is possible to constrain the 
mass (m g ) of this particle. 

In section 2 we present the new solution for the Yukawian gravitational potential 
for a thin disk and compare it, in section 3, with the Newtonian solution obtained by 
[5]. In section 4 we present the main results and conclusions. 



2 The rotation curve for a thin exponential disk for a Yukawian 
gravitational potential 

Toomre [5] showed that the potential of a thin disk can be written as follows 

roo 

(j>{r,z) = 2-kG J (kr)S(k)e~ klzl dk, (5) 
JO 

where S(k) is related to the surface density, n(r), through the following Bessel integral 

K r ) = / J (kr)S(k)kdk, (6) 
Jo 

where 

/■OO 

S(k) = / J (ku)fi(u)udu (7) 



comes from the Fourier-Bessel integral theorem. 
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It is straightforward to show that Eq. ([5]) satisfies the Newtonian field Equation. 
Once fi(r) is given one can obtain the potential of the disk, and from the centrifugal- 
equilibrium condition, namely, 



z=0 



one obtains the rotation curve of the disk. 

Freeman ([5]) applied the above equations to obtain the rotation curve for an 
exponential disk, whose surface density, n(r), in cylindrical coordinates is given by 

p(r) = ^e- r/ \ (9) 

where b is the scale length of the disk. 

It is worth stressing that the motivation for adopting an exponential disk comes 
from observations. Surface photometry shows that the two main components of most 
spiral and SO galaxies are spheroidal and disk components. 

In particular, the radial surface brightness distribution of the disk follows an 
exponential law, namely, I(r) — Ioe~ r ^ b . This implies that the surface density is 
exponential too, as we are considering. 

Finally, the circular velocity, i.e., the rotation curve, for this exponential disk in 
the Newtonian gravitation reads 

v 2 (r) = 7vG/i b (J-) 2 (I K - hKi\ (10) 

where I and K are the modified Bessel functions, which are calculated at r/2b. 

We now consider how to get a solution for the rotation curve of a thin disk for a 
Yukawian potential, using the Toomre approach. 

The Yukawian potential for a thin disk can be obtained by rewriting Eq. (JS| 
properly, namely, 



2-kG 



n 



Z (\/\k 2 - A- 2 r)S(k)er k ^dk, (II) 



where 



W ' ' ' JnMfc 2 -A-2| r ) 



-1 

(Vlfe 2 - A~ 2 | r) for k < A -1 



(12) 



As in the Newtonian case the Bessel functions appear in the solution presented 
above. To account for the extra term u (f>/\ 2 " in the Poisson equation the argument of 
the Bessel functions must now contain "A". 

Obviously, the direct substitution of the above equations satisfy the Yukawian 
gravitational potential Eq. (|4)l. 

We now apply the above equations to obtain the rotation curve for an exponential 
disk, whose surface density, n(r), in cylindrical coordinates is given by Eq. (|9]). 
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After some manipulation and an eventual change of variables the rotation curve 
reads 

v 2 (r) — 

^/x 2 -b 2 /X 2 



Jl ( T -^x?- -6 2 /A 2 ) dx 



Jb/X (l+*2)3/2 

(1 + 3,2)3/2 l ybV ' ) 



(13) 



It is worth stressing that we have parameterized the above equation in terms of the ratio 
"b/X" , which helps one to analyze how different a theory with a Yukawian gravitational 
potential is as compared to the Newtonian theory. 

Note that for b/X <C 1 the second integral of the Eq. (I13p goes to zero and the first 
one becomes 



*Ji(tF) 
o (1 + z 2 ) 3 / 2 



dx - 



1 r 

2 b 



{IqKq-IiKx). 



(14) 



which substituting into Eq. (|13f) gives, as expected, the Newtonian rotation curve. 

In the next section we integrate numerically the Yukawian rotation curve for 
different values of "6/ A" and compare them with the Newtonian rotation curve. 



3 Yukawian versus Newtonian disk 

We now compare the Newtonian rotation curve with those that come from a Yukawian 
gravitational potential. To proceed, however, one needs to integrate Eq. (|13p . which 
has three parameters to be specified, namely, A, [iq and b. 

The A parameter is in principle free, since we are not considering any specific 
Yukawian gravitational theory. The parameters fiQ and 6 can be inferred from 
observations, as already mentioned. For reasons that become clear later on we here 
only need to know b. 

There is in the literature a series of studies concerning the determination of the 
disk parameters of spiral and SO galaxies. It is worth noting that there is not any 
universal value to the scale length b, different galaxies present different values for this 
parameter. From these studies one can conclude that b ~ 1 — 10 kpc [7]. 

Since A is a free parameter, whose value depends on the particular theory adopted, 
we do not fix any particular value for it. Instead, since the scale length is constrained 
by observations, we present the rotation curves parameterized in terms of "b/X". 

In Fig.[T]we compare the Newtonian rotation curve with different Yukawian rotation 
curves. We plot the velocity, in terms of y/G fio b, which has dimension of velocity, versus 
the distance to the center of the disk in units of the scale length, 6. 

The Newtonian rotation curve has a maximum velocity around r/b ~ 2, and for 
larger values of r the velocity decreases monotonically. This is a well known result. 
On the other hand, the behavior of the Yukawian rotation curves is the following: the 
velocity increases up to a maximum value and then decreases; the greater the ratio 
"b/X" , the lower the maximum velocity is. The size of the disk is strongly dependent 
on the value of the ratio "b/X" . The greater the ratio "b/X" , the smaller the disk is. 
As expected, for b/X <C 1 the rotation curve becomes Newtonian. 
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Fig. 1 Comparison between the Newtonian rotation curve with different Yukawian rotation 
curves parameterized in terms of "6/ A". 



These behaviors of the Yukawian rotation curves strongly constrain the value of A. 
As a result, for theories where A is the Compton wavelength, strong constrains can be 
imposed on the putative mass of the graviton. 

To the Yukawian rotation curve be approximately the Newtonian one, we need 
b/X <C 1. Recall that the Compton wavelength, A c , of a particle of mass m g reads 

A - h 

Ac — , 

m g c 

where c is the velocity of light and h is the Planck constant. For b — 10 kpc we should 
have 

m g « l(T 59 g. 
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4 Discussions and conclusions 

We present a new solution for the rotation curves of a thin disk for a Yukawian 
gravitational potential. This thin disk can represent the disks present in spiral and 
SO galaxies. It is important to stress that the mathematical tool developed here can 
be used to study rotation curves in any model or theories of gravity which produce a 
gravitational potential of the Yukawa type. 

In particular, we compare the Yukawian rotation curve with the Newtonian one. 
The Newtonian rotation curve has a maximum velocity and for larger values of r the 
velocity decreases monotonically. 

The main characteristics of the Yukawian rotation curves are the following: the 
velocity increases to a maximum value and then decreases; the greater the ratio "b/X" , 
the lower the maximum velocity is. The greater the ratio "6/A", the smaller the disk 
is. As expected for b/X <C 1 the rotation curve becomes Newtonian. 

These behaviors of the Yukawian rotation curves strongly constrain the value of A. 
As a result, for theories where A is the Compton wavelength of a exchange particle, say 
a graviton, strong constraints can be imposed on its putative mass. For 6 = 10 kpc, for 
example, we find m g <C 10 _59 g. 

It is worth stressing that the best bound on the graviton mass from planetary 
motion surveys is obtained by using Kepler's third law to compare the orbits of Earth 
and Mars, yielding m g < 10 -54 g [8]. With such a mass and for b = 10 kpc one would 
obtain that b/X ~ 10 5 , which would imply that galactic discs could not exit (see Fig.[T]) 
for a Yukawian gravitational potential with this value of m g . 

A study of interest, concerning a Yukawian gravitational potential, would be 
focused on a complete model for the rotation curve including, besides a disk, a halo. 
Such a study might impose additional constraints on A or on the putative mass of a 
graviton. 

Also interesting is to study clusters of galaxies in Yukawian gravitational theory. 
Since the size of clusters extends to ~ 1 Mpc, strong limits may be imposed on A or on 
the mass of the graviton. We leave these studies, however, to another papers to appear 
elsewhere. 

Although we leave these studies to another papers, our results summarized by the 
a parameter in Fig. 1, can already be compared to the galaxy rotation curves of the 
sample studied by, for example, Ref. [5]. 

In general, the behavior of the rotation curves, which come from the observations, 
is the following: the velocity increases as a function of the distance r to the center of 
the galaxy; it then becomes nearly flat for large values of r. 

The innermost part of the rotation curves, where the velocity is an increasing 
function of r, is dominated by the galactic disk. This very behavior allows us to compare 
our results with observations. On the other hand, the outermost part of the rotation 
curve, where it becomes nearly flat, is dominate by the galactic halo (see, e.g., [10] , 
page 371). 

In particular, Brownstein and Moffat [9] analyzed the rotation curves of a sample of 
101 galaxies compiled from [TT„ 12 13. 14 15 16 . All these galaxies present a behavior, 
in the innermost part of their rotation curves, which are consistent with a Newtonian 
gravitational potential. There is, as a result, a perfect agreement between our results 
and the observations if a <C 1. This very fact, as already mentioned, strongly constrain 
the value of the graviton mass. 
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Moreover, numerous studies of the Tully-Fisher (TF) relations and their 
applications have been conducted in the past. In most cases their aim was to find 
observables that reduce the scattering and thus improve this relation as a tool for 
measuring distances to spiral galaxies, and also to gain insights in the formation and 
structure of galaxies [17] . 

In particular, il4] investigate the statistical properties of the Tully-Fisher (TF) 
relations for a volume-limited complete sample of spiral galaxies in the nearby Ursa 
Major Cluster. They concluded that the TF relation reflects a fundamental correlation 
between the mass of the halo and the total baryonic mass of the galaxies. This kind of 
result certainly reinforces our interest to apply the formalism here discussed on a more 
complete model for rotation curves which includes a massive dark halo. 

Acknowledgements JCNA would like to thank the Brazilian agency CNPq for partial 
support. 



References 

1. Moffat, J.W., Sokolov, I.Yu.: Phys. Lett. B 378, 59 (1996) 

2. Piazza, F., Marinoni, C: Phys. Rev. Lett. 91, 141301 (2003) 

3. Rodriguez-Meza, M.A., et al.: Gen. Relativ. Gravit. 37, 823 (2005) 

4. Signorc, R.L.: Mon. Not. R. Astron. Soc. 364, 1219 (2005) 

5. Freeman, K.C.: Astrophys. J. 160, 811 (1970) 

6. Toomre, A.: Astrophys. J. 138, 385 (1963) 

7. Flores, R., Primack, J.R., Blumenthal, G.R., Faber, S.M.: Astrophys. J. 412, 443 (1993) 

8. Larson, S.L., Hiscock, W.A.: Phys. Rev. D, 61, 104008 (2000) 

9. Brownstein, J.R., Moffat, J.W.: Astrophys. J. 636, 721 (2006) 

10. Peacock, J. A.: Cosmological Physics. Cambridge University Press, Cambridge (1999), p. 
371 

11. Begeman, K., Broeils, A.H., Sanders, R.H.: Mon. Not. R. Astron. Soc. 249, 523 (1991) 

12. Sander, R.H.: Astrophys. J. 473, 117 (1996) 

13. de Blok, W.J.G, McGaugh, S.S.: Astrophys. J. 508, 132 (1998) 

14. Verheijen, M.A.W., Sancisi, R.: Astron. Astrophys. 370, 765 (2001) 

15. Some, Y.: Astrophys. J. 458, 120 (1996) 

16. Romanowsky, A. J., et al.: in IAU Symp. 220, Dark Matter in Galaxies, ed. S. Ryder, et 
al. (San Francisco: ASP), 165 (2004) 

17. Hinz, J.L., Rix,H.-W., Bernstein, G.M.: Astrophys. J. 121, 683 (2001) 



